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Abstract. For a d-dimensional polyhedral complex P, the dimension of the 
space of piecewise polynomial functions (splines) on P of smoothness r and 
degree k is given, for k sufficiently large, by a polynomial f[P, r, k) of degree 
d. When d = 2 and P is simplicial, in [T] Alfeld and Schumaker give a formula 
for all three coefficients of /. However, in the polyhedral case, no formula is 
known. Using localization techniques and specialized dual graphs associated to 
codimension— 2 linear spaces, we obtain the first three coefficients of f(P, r, k), 
giving a complete answer when d = 2. 



1. Introduction 

In [5] , Billera used methods of homological and commutative algebra to solve a 
conjecture of Strang on the dimension of the space of C 1 splines on a generic planar 
simplicial complex. The algebraic approach to the study of splines was further 
developed in work of Billera and Rose [4] , [5] , and Schenck and Stillman [16] , [17] . 
For a (not necessarily generic) triangulation A of a simply connected polygonal 
domain having /{* interior edges and /q interior vertices, Alfeld and Schumaker 
showed that for k > 3r + 1, the dimension of the vector space C£(A) of splines on 
A having smoothness r and degree at most k is given by: 

Theorem 1.1 ([2]). 

where a = ^2 a i> a i = max{(r + 1 + j(l — n(vi))), 0}, and n{vi) is the number 
of distinct slopes at an interior vertex Vi . 

In higher dimensions, a spectral sequence argument [13j shows that the previous 
theorem continues to hold, in the sense that for k ^> 0, the dimension of C£(A) is 
given by a polynomial in k (the Hilbert polynomial) , and the coefficients of the three 
largest terms of this polynomial are determined by a higher dimensional analog of 
the formula given by Alfeld and Schumaker in the planar case. 

In this note, we study splines on a polyhedral complex P. In general, polyhedral 
complexes do not lend themselves to the full range of techniques available in the 
simplicial case, and so have been studied less. In [T5], Schumaker obtains upper 
and lower bounds in the planar case, and in 0], Billera and Rose obtained the first 
two coefficients of the Hilbert polynomial. It is possible to consider all the C r k (P) 
at once by passing to a module over a polynomial ring, and in [20] Yuzvinsky uses 
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sheaves on posets to obtain results on the freeness of this module. In [IT], [12], Rose 
studies cycles on the dual graph of P. Our key technical innovation is a refined 
version of the dual graph, depending on a choice of codimension-2 linear space; 
used in conjunction with localization techniques. 

The main result of this paper is a formula for the third coefficient of the Hilbert 
polynomial, in the case of a polyhedral complex. In particular, for a planar polyhe- 
dral complex, our result gives a formula for the dimension of CJ1(P) for all k ^> 0. 

Example 1.2. Let P be the polygonal complex depicted below: 




The following table gives the dimensions for CT(P) for small values of r, assuming 
k 3> 0. For a planar P, Corollary 13. 141 gives a formula for f(P,r,k) depending on 
the combinatorics and geometry of P, "explaining" the values below. 



r 


dim* C£(P) 





2k 2 + 2 


1 


2k' 2 -6fc + 10 


2 


2k 2 - 12fc + 32 


3 


2k 2 - 18fc + 64 


4 


2k 2 - 24k + 110 



Throughout the paper, our basic references are de Boor [7] (for splines) and 
Eisenbud [8] (for commutative algebra). 

2. Algebraic preliminaries 

Definition 2.1. The dual graph Gp of a d-dimensional polyhedral complex P is a 
graph whose vertices correspond to the d-faces of P, with an edge connecting two 
vertices iff the corresponding two d-faces meet in a (d — l)-face. 

Definition 2.2. The star of a face a of P is the set of all faces t £ P such that 
there exists a face v € P containing both r and a. 

Put another way, the star of a is the smallest subcomplex of P containing all 
faces which contain a. Let P be a <i-dimensional polyhedral complex embedded 
in M. d , such that P is hereditary; this means for every nonempty face a of P, 
the dual graph of the star of a is connected. Basically, P should be visualized 
as a polyhedral subdivision of a manifold with boundary. A C r -spline on P is 
a piecewise polynomial function (a polynomial is assigned to each d-dimensional 
cell of P), such that two polynomials supported on d-faces which share a common 
(d— l)-face r meet with order of smoothness r along that face. We use P? to denote 
the set of interior i faces of P (all d-dimensional faces are considered interior), and 
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let f° = \Pf\. The set of splines of degree at most k (i.e. each individual polynomial 
is of degree at most k) is a vector space, which we will denote C£(P). 

2.1. Splines and syzygies. The first important fact is that the smoothness con- 
dition is local: for two d-cells o~\ and &i sharing a common (d— l)-face r, let l T be a 
nonzero linear form vanishing on r. Billera and Rose ([4], Corollary 1.3) show that 
a pair of polynomials supported on <jj,i = 1,2 meet with order r smoothness 
along r iff 

- /a- 

In [IS] , Schumaker gave a formula for the dimension of C£(P) in the planar case, 
when P is the star of a vertex, as depicted in Example 12.31 (note that such a P will 
always be simplicial). Schumaker observed that a spline is a syzygy on the ideal 
generated by the r + 1 st powers of linear forms vanishing on the interior edges. To 
see this, consider the example below: 

Example 2.3. A planar P which is the star of a single interior vertex vq at the 
origin. 



(-LD 



(0,1) 




(1,0) 



(-1.-D 

Beginning with the simplex in the first quadrant and moving clockwise, label the 



■+1 



polynomials on the triangles /i , . . . 

aiy 

a3(x + yy 
a 4 x r+1 



To obtain a global C r function, we require 
= h-h 

= h — h 
1 = 73-/4 

= h-h 



r+1 _ 



0, a syzygy on the 



Summing each side above yields the equation 
ideal 

I=(y r+1 ,(x-y) r+1 ) (x + y) r+1 ,x r + 1 ). 
It is easy to see that the process can be reversed. The upshot is that if P is the 
star of a vertex, then the spline module consists of K[x, y] syz(I); the summand 
y] corresponds to the "trivial" splines /1 = f% = h = h- Hence, splines are 
intimately connected to commutative algebra. 

2.2. Graded modules, Hilbert polynomial and series. As first observed by 
Billera and Rose, one way to study the dimension of the vector space C^(P) is to 
embed P in the hyperplane Xd+i = 1 C R d+1 , and form the cone P over P, with 
vertex at the origin. If C^(P) is the set of splines on P such that each polynomial 
is homogeneous of degree k, then Billera and Rose show ([4], Theorem 2.6) there is 
a vector space isomorphism between C^(P) and C^(P). It is easy to see that the 
set of splines of all degrees 

c r (P) = |J cx(p) 

fc>0 



4 



TERRY MCDONALD AND HAL SCHENCK 



is a graded module over the polynomial ring R in d + 1 variables. It follows from 
basic commutative algebra (see [14]) that the dimension of all the graded pieces of 
C r (P) is encoded by the Hilbert series: 

HS(C r (P),t) = £dimC£(P)t fe = (i 9 ® )d+1 - 

for some g(t) £ Z[t], and that for fc ^> 0, the dimension of C£(P) is given by the 
Hilbert polynomial HP(C r (P), k), which is an element of Q[k] of degree d. An easy 
induction shows that 

HS(R, t) = and HP(R, k) = + * 

We need one last bit of notation: R(—i) will denote the polynomial ring R, consid- 
ered as a graded module over R, but with generator in degree i: 

HS{R{-i),t) = {1 _ f t)d+1 and HP(R(-i), k) = ^ + J " 1 

Lemma 2.4. [Billera-Rose, [4]] Let P be ad-dimensional polyhedral complex. Then 
there is a graded exact sequence: 

— ► C r (P) — > R fd © R f *-i(-r - 1) R^- 1 — > N — > 

11 



where <f> = 



jr-\ 

7"; j 



Write [<?d | -D] for <\>. To describe note that the rows of da are indexed by 
r G P*d-\- If cr i' CT 2 denote the d— faces adjacent to r, then in the row corresponding 
to r the smoothness condition means that the only nonzero entries occur in the 
columns corresponding to <7x,0"2, and are ±(+1, —1). When P is simplicial, dd is 
the top boundary map in the (relative) chain complex. 

A main result of [4] is that N is supported on primes of codimension at least 
two. Combining this with the exact sequence of Lemma 12 . 41 allows Billera and Rose 
to determine the first two coefficients of the Hilbert polynomial (though the result 
is phrased in terms of the Hilbert series). 

3. Main Theorem 

We begin by sketching our strategy. It follows from additivity of the Hilbert 
polynomial on exact sequences and Lemma 12.41 that obtaining the coefficient of 
k d ~ 2 in the Hilbert polynomial of C r {P) is equivalent to obtaining the coefficient 
of k d ~ 2 in the Hilbert polynomial of N. Since 

iV^(0 R/l r T +1 )/d d , 

tG P° 

every element of N is torsion. Using localization, we first show that the 
codimension-2 associated primes of N must be linear, then give a precise descrip- 
tion of which codimension-2 linear primes actually occur. This leads to an explicit 
description of the submodule of N supported in codimension-2. Elements of this 
submodule are the only elements of N which contribute to the k d ~ 2 coefficient of 
the Hilbert polynomial, and the formula follows. 
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3.1. The codimension two associated primes of N. 

Lemma 3.1. Any codimension-2 prime ideal J associated to N contains a linear 
form l T , for some r G P^—i- 

Proof. From the description 

iv^(0 R/r+i)/a d , 

tPP° 

it follows that if no l T is in J, then all the l T are invertible in Rj, so that Nj 
vanishes. □ 

Lemma 3.2. Let £ be a codimension-2 linear space. If a 6 Pd has at most one 
facet whose linear span contains then every generator of N corresponding to a 
facet of a is mapped to zero in the localization . 

Proof. In Ri(£), any l T such that £ % V(l T ) becomes invertible. As N is the cokernel 
of 

rr+l 

d d 

ir + l 

in the right hand submatrix D of 4>, all the forms (or all save one) V T +1 such that r is 
a facet of a become units. As the column of the left hand (dd) matrix corresponding 
to cr has nonzero entries only in rows corresponding to facets of er, this means that 
every generator corresponding to a facet of a has zero image in the localization, 
and the result follows. □ 

Theorem 3.3. Any codimension-2 prime ideal J associated to N is of the form 
(l Tll l T2 ) for Ti G Pj_i such that V(l Tl ,l T2 ) has codimension-2. 

Proof. If there do not exist two k as above, then by Lemma l3~Tl V(J) is contained 
in exactly one hyperplane which is the linear span of t G P®-i- Thus, in Rj, all but 
one of the l T become units, and the proof of Lemma l3~2l shows that Nj vanishes. □ 

Theorem [33] gives an explicit set of candidates for the codimension two primes of 
AT. As noted earlier, in 4 , Billera and Rose showed that all associated primes of N 
have codimension at least two (this also follows from the arguments above), so the 
theorem identifies all candidates for the associated primes of minimal codimension. 

Proposition 3.4. Let N be a finitely generated, graded R = K[xi, . . . , xj] module 
of dimension m < d. Let V be the set of minimal (codimension d — m) associated 
primes of N, and for Q G V , define 

N(Q) = {n G N | ann(n) C Q}. 
IfHP(N,k) = a m k m + 0(k m - 1 ), then 

HP(® N(Q),k)=a m k m + 0(k m - 1 ). 

Proof. First, note that N(Q) N. Then the direct sum of the N(Q), taken over 
the finite set of Q G "P, maps to N, and so there is an exact sequence 

— v K — > N(Q) — > N — > C — > 0. 
Qev 
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By construction, C is supported in codimension greater than d — to, so has 
Hilbert polynomial of degree at most m — 1. Next we show that N(P)q — for 
any P,Q e V such that P ^ Q. From the definition of N(P), for any n <E N(P), 
J ann(rt) = P. As Q ^ P, there exists p E P,p ^ Q. Since p m ■ n = for 
some m, and since p is a unit in Rq, we see that N(P)q = 0. Localization is an 
exact functor, and so localizing the exact sequence above at Q E V yields an exact 
sequence: 

— K Q — 7V(Q)q — > iVq — C Q — > 0. 
Localizing the exact sequence — ► N(Q) — ► iV yields an exact sequence 

— ► N(Q) Q — > TVq, 

which implies that i^Q must vanish. Thus, -ftT is also supported in codimension 
greater than d — to, so has Hilbert polynomial of degree at most m — 1. The result 
follows by additivity of Hilbert polynomials in an exact sequence. □ 

3.2. The graph associated to P and a codimension two subspace. In this 
section, we analyze the codimension-2 associated primes of N in greater depth. 
In the simplicial, planar case, the Alfeld-Schumaker formula shows that there is a 
contribution to the constant term of the Hilbert polynomial from the local behaviour 
at interior vertices. As noted in Example 1 2. 3 [ a syzygy on the linear forms adjacent 
to a single fixed vertex corresponds to a loop around that vertex. 

Example 3.5. For the planar polyhedral complex P appearing in Example II. 2 [ let 
Vi be the vertex associated to the face labelled fi. The dual graph G{P) is: 




In [T2], Rose shows that if the dual graph Gp of P has a basis of disjoint cycles, then 
the projective dimension and Hilbert series of C r (P) are determined by the case 
when Gp has a single cycle, which Rose analyzed in [11]. We now define a refined 
version of dual graph, which depends on P and the choice of a codimension-2 linear 
subspace. 

Definition 3.6. Let P be a d- dimensional polyhedral complex embedded in R d ; 
and £ a codimension-2 linear subspace. G^ (P) is a graph whose vertices correspond 
to those a S Pd such there exists a (d — l)-face of a whose linear span contains 
£. Two vertices of G$(P) are joined iff the corresponding d-faces share a common 
(d — l)-/ace whose linear span contains £. 



PIECEWISE POLYNOMIALS ON POLYHEDRAL COMPLEXES 



7 



Example 3.7. We continue to analyze the complex P of Example 11.21 For each 
interior vertex v of P, G V {P) consists of a triangle. However, there is another 
codimension-2 space to consider: if £ is the point at which (the linear hulls of) the 
three edges connecting interior vertices to boundary vertices meet, then G^{P) is 




Let P' be obtained by moving the top vertex of P a bit to the right. Then P 
and P' are combinatorially equivalent, but in P' there are no sets of > 3 concurrent 
V(I t ),t S PS-i> except at the interior vertices. In particular, G^(P') is acyclic. 

Lemma 3.8. For any a £ Pd, there are at most two facets of a whose linear spans 
contain a given codimension-2 linear space £. 

Proof. Suppose the linear spans of three facets T\ , T2, T3 of a meet in a codimension- 
2 linear space £. For each V(l n ), a lies on one side of the hyperplane; so a lies 
between V(l Tl ) and V(l T2 ). Since V(l T3 ) contains £ = V(l Tl ) H V(l T2 ), this means 
V(l T3 ) would split a, a contradiction. □ 

Corollary 3.9. G^(P) is homotopic to a disjoint union of circles and segments. 

Proof. By Lemma [5T51 the valence of any vertex v € G^(P) is at most two. □ 

Theorem 3.10. For a polyhedral complex P and codimension-2 linear prime 

N m)^ (R/i*)m 

where ip £ Hi(G^(P)) means %j} is a component of G$(P) homotopic to S , and 

Iij, = {l r T +1 I t S Pd-\ corresponds to an edge oftjj). 

Proof. By Corollary 13. 91 G$(P) consists of a disjoint union of cycles and segments. 
By Lemma 13.21 all generators of N which lie in a segment are mapped to zero in 
the localization Nj^y For each d— face a corresponding to a vertex in a cycle, 
note that there are two (d — 1)— faces t\, ti of a such that l Tl , l T2 are not units in 
Rite)] every other linear form defining a facet of a becomes a unit. Reducing the 
column of dd corresponding to a by the columns of -D/(£) having a unit entry gives 
a column with nonzero entries only in rows corresponding to t\ and T2 . Repeating 
the process shows that the cycle corresponds to a principal submodule of iV/(£), 
with the generator quotiented by the [r + l) s * powers of the forms corresponding 
to the edges of the cycle. □ 

Corollary 3.11. Let Q be a codimension-2 associated prime of N with V(Q) = £, 
and let N(Q) C N be as in Proposition \3.4\ Then 

JV(Q)~ R/I4, 

VGffi(G 5 (P)) 
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where the direct sum on the right hand side is as in Theorem \3.1(A 

Proof. Let K(£) be the submodule of N generated by r e P d -% such that £ is 
not contained in the linear span of r. These generators correspond exactly to the 
columns of D in which l T is a unit in the localization at Q. Localize the exact 
sequence 

— ► K{C) — > N — > N/K{C) — ► 0. 
at Q. By construction, 

(N/K(£)) q ~N(Q) q ~(@N(P))q. 

So we have that N(Q) ~ N/K(£). To finish, we need to show that N/K(£) de- 
composes as the direct sum (over cycles in G^(P)) of the R/I^. In N/K(£), we 
may use the columns in D which contain a unit to reduce each column of d d - After 
reducing, a column corresponding to a G P c i has at most two nonzero entries. If 
a is part of a segment in G^(P), then in N/K(£) the generators corresponding to 
t C a reduce to zero. Otherwise, a corresponds to a vertex in a cycle of G$(P), and 
has exactly two (modulo K(£)) nonzero entries. So in each cycle, all the generators 
are equivalent in N/K(£), and are quotiented by the appropriate i£ + . □ 

Theorem 3.12. Let M be a graded module, and let ad- 2 {M) denote the coefficient 
of k d ~ 2 in HP(M, k) . As in Proposition \3.4\ let V be the set of codimension-1 
primes associated to N, and for Qi £ V let £j = V(Qi). Then ad- 2 (C r (P)) = 

a d - 2 {R U - f *- 1 ) + a d - 2 {R{-r-l) f ^)+ Yl a^R/I^). 

Proof. By Lemma l2.4l and the additivity of Hilbert polynomials on exact sequences, 

a d - 2 (C r (P)) = a d -2(R fd - f «- 1 ) + a d . 2 (R(-r - 1)^) + a d - 2 {N). 
It follows from Proposition OQ] that 

a d - 2 {N) = J2 a d -2(N{Qi)). 

QiEV 

By Corollary GUU 

JV(Qi)~ R/I^, 
^•eKi(G ?i (p)) 

and the result follows. □ 

The value of a d - 2 (R(— m)) is a Stirling number of the first kind, and can be 
written out explicitly. In order for Theorem 13.121 to be useful, we need to know 
a-d- 2 (R/Iti>), which is provided by the following lemma: 

Lemma 3.13. Let 1^ = (t^~ , . . . , C K[xo, ■ ■ ■ , x d ] be a codimension-2 ideal, 
minimally generated by the n given elements. Define 

atyO=LSlJ> s 1 (rJj) = (n-l)a(Tp)+n-r-2, s 2 {^) = r + l-{n-l)a{^). 

Then the minimal free resolution ofR/I^ is: 

R(-r-l-a(tp)) Sl ^ 
— ► © — >R(-r-l) n — >R — >R/I i> — >0. 

R{-r-2-a(ip)) S2 ^ 
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Proof. See Theorem 3.1 of |16j : the key step involves showing that a certain matrix 
has full rank, which was established by Schumaker in 1SJ . □ 

It follows from Lemma T3. 131 that the Hilbert polynomial of R/I^ is given by: 
k+d\ fk+d-r-l\ , ^(k+d-r-l-a.(ib)\ , . fk+d-r-2~a{ib) 



Corollary 3.14. If P is a hereditary planar polyhedral complex, then 
HP(Cr(P), k) = lie + 3/2-2(r + l)/f fc + + yj + £ 



V> 3 '£i?i(G 4i (P)) 



where 



2 



2r + 3 + a(^) - n(l + a(^-))) ■ 



4. Examples and connection to simplicial case 

We close with some examples, and a discussion of the relation to the simplicial 
case. We begin by applying Corollarv l3.14l to Example 1 1.2 1 As we saw in Example 
13. 7\ there are four £ at which Hi(G^(P)) ^ 0, and each 1^ has three generators. 
Hence the Cj are all the same, and equal to 

( r 2 2 ) + ^ (2r + 3 + a(^) 3(1 + afcfc))) . 

which simplifies to 

fr + 2\ ,r + l./ ,r + l.\ 

( 2 J + L— J(r-L— J> 

Applying Corollary 13 . 141 yields : 



r 


dim K C r k {P) 


h h 2 , 3/ 2 -2(r+l)/; J , 
2 2 K 


h + ((l)-i)f? 


4((''+ 2 )+a(r-a)) 





2k' 2 + 2 


2fc 2 


-2 


4 


1 


2k' 2 -6fc + 10 


2fc 2 - 6fc 


-2 


12 


2 


2fc 2 - 12/s + 32 


2fc 2 - 12k 


4 


28 


3 


2fc 2 - 18/fc + 64 


2k' 2 - 18k 


16 


48 


4 


2fc 2 - 24/c+ 110 


2k' 2 - 24fc 


34 


76 



As in Example 13. 7\ consider the configuration P' obtained by perturbing a vertex 
in Example 11.21 so the three edges defining £ no longer meet. Then there are only 
three nontrivial Cj , and we have 
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r 


diniffi CX(P') 


h h 2 | 3/ 3 -2(r+l)/ 1 u , 


/ 2 +(©-i)A° 


3(ff)+a(r-a)) 





2k' 2 + 1 


2fc 2 


-2 


3 


1 


2k' 2 -6k + 7 


2fc 2 - 6k 


-2 


9 


2 


2k' 2 - 12k + 25 


2k' 2 - 12k 


4 


21 


3 


2fc 2 - 18/c + 52 


2k' 2 - 18fc 


16 


36 


4 


2fc 2 - 24/c + 91 


2k' 2 - 24k 


34 


57 



Example 4.1. In this example, we show that G^(P) can have several disjoint 
cycles. Let P be as below: 




If £ corresponds to the central vertex, then G^{P) is: 





For these two cycles, 1^ has only two generators, and the Cj value for such an ideal 
is always (r + l) 2 . There are four cycles for which Jy, has three generators, so by 



Corollary 13 . 141 we have (this table omits column for &k 



3/2-2(r+l)/° 



k): 



r 


dim R q;(P) 


/ 2 +(©-i)A° 


4{( r f)+a{r-a)) 


2(r+l) 2 





4fc 2 + 2 


-4 


4 


2 


1 


4k' 2 - 12fc+16 


-4 


12 


8 


2 


4k' 2 - 24k + 54 


8 


28 


18 


3 


4k' 2 - 36fc + 112 


32 


48 


32 


4 


4fc 2 -48fc+194 


68 


76 


50 
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Example 4.2. In our final example, we look at a honeycomb configuration, that 
is, the polyhedral complex consisting of seven hexagons, pictured below: 



For each interior vertex v, G V (P) will be a triangle and a line segment; while if 
£ is the point at the center of the diagram, Gt(P) is a hexagon. However, each of 
the lip is generated by three elements, and so Corollary 13. 141 yields: 



r 


dim R C r k {P) 


/2L2 , 3/ 2 -2(r+l)/;' , 
2 2 K 


/ 2 +(©-!)A° 


7(( r f)+a(r-a)) 





7fk Tfk ~\~ 2 




-5 


7 


1 






-5 


21 


2 


~2 k ~2~ k ~\~ 56 


7 bl r. 

2 2 


7 


49 


3 


77 k — ~~y~k ~\~ 115 


V>,2 VbjL 


31 


84 


4 


JJfc a -4 fc + 200 


7 u'2 #J 

2 K 2 K 


67 


133 



4.1. Connection to the simplicial case. In the case where P is actually a sim- 
plicial complex, the first point to notice is that when two facets n , Ti of a G Pd 
meet, then they actually meet in a face of a, necessarily of codimension-2. Thus, 
in the simplicial case, the only £ such that Hi(G^(P)) ^ are £ € P®-2- For a 
planar P which is simplicial, there are /q such faces. The formula for the value Cj 
which appears in Corollary 13. 141 can be rewritten as: 

+ 2^ 



where <x; is as in the Alfeld-Schumaker formula. So in the planar, simplicial case, 
the formula of Corollary 13. 141 is the polynomial appearing in Theorem 1.1. 

Concluding Remarks and Questions: 

(1) To obtain the next coefficient of the Hilbert polynomial will be difficult; 
indeed, even in the simplicial case no general formula is known. 

(2) As noted in [T5], the value for which the dimension of C£(P) becomes 
polynomial is the Castelnuovo-Mumford regularity the sheaf associated to 
C r (P); it would be interesting to determine the regularity of C r (P). 

(3) It is possible to generalize these results to the case of mixed smoothness, 
which was considered for planar simplicial complexes in [S]; we leave this 
to the interested reader. 

(4) In the simplicial, planar case, [6] shows that C r (P) free — > C r ~ 1 (P) free. 
Is this true for P polyhedral? 

Acknowledgments: Macaulay2 [TU] computations were essential to our work. 
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